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Hans- Jurgen Kayser publications 

published by bk-teachware - “the“ CAS-Specialist 

^ Elektronische Arbeitsblatter zur Stochastik in der Sek I mit Derive, SR-27 
Derive im Mathematikunterricht der Sek 1, SR-35 
^ Analysis mit Derive, Lehr- und Arbeitsbuch fur die Jahrgangsstufe 11, SR-38 
^ Lineare Algebra und Geometrie mit Derive, SR-42 

Derive im Analysisunterricht der Jahrgangsstufen 12 und 13, SR-45 
Derive im Stochastikunterricht der Sekk II, SR-45 

Visit http : / / shop . bk-teachware . com 

A Website with a bundle of Derive-materials from Maria Koth, Vienna: 

http : //www.mathe-online . at/lernpfade/derive/ 


GREAT NEWS from Virgin Islands 
2 new Calculator APPS are now available from the author of 

Calculus Made Easy: 

1) TRIGONOMETRY MADE EASY - 

Solves ANY Trigonometry Problem instantly. 

Solve any triangle completely, review all trig formulas, graph and model with sin(x)-functions, 
Unit Circle coordinates, Convert degree/radian, revolutions/degrees and much more. 

2) MATRICES AND VECTORS MADE EASY - 
Solves ANY Matrix and Vector problem instantly. 

Edit and View Matrices, do 20+ matrix computations such A+B, A*B, Inverse, Determinant, 
Eigenvalues, 

LU-Factorization, Solve systems of equations by viewing steps, Apply Matrices. Add and Graph 
Vectors, 

Compute Dot/Cross Products, UnitVectors, and much more. 


Calculus Made Easy 6.0 is now available. 

Major Upgrades make CME now even more powerful and easier to use: 

1) NEW: Perform Logarithmic Differentiation. Find Average Rate of Change. Simpson Rule 
for 

area approximation. 

2) NEW: Do both Fundamental Theorems of Calculus. 

3) NEW: Graph ...functions, ...Average Value, ...Area between Functions, Under Curves, 
two Areas of Equal Size, Net and Total Area. 

4) NEW: Display setup of Integrals before actual Integration. More Information on Differential 
Equations: Logistic Growth, Exponential Growth, etc. 

5) NEW: Module to find Extrema of 2-dim. functions significantly improved. 

6) NEW: Take a break from Calculus and play a game of Yahtzee. 

7) And much more...i.e. scroll answers too wide for screen. 

Available at www.ti89.com . Trial Versions available. 
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Dear DUG-Members, 

First of all I'd like to apologize once more for the 
delay of publishing this DNL#60. As I wrote in my 
last information mail we were in Namibia for three 
weeks in December. But this was not a mathless 
time for me because when my travel companions 
had heard that I was a math teacher they and I dis- 
covered a lot of mathematics during the trip. So I 
remembered David Halprin's great “Sand Dunes” 
article climbing up the dune in the Namib desert, 



quiver trees reminded me on fractal structures, 




a window of the church in Windhoek invites to be 
modeled by geometric figures 



and even the horns of a rhino let us think on 
mathematical curves. 



I intended to include Peter Schofield's and Don 
Phillips' articles but as you can see this DNL is 
overfull again. I am very proud that Bernhard 
Kutzler delivered a lecture in written form which he 
gave last summer in Atlanta and that Josef Lechner 
presents his rose in winter time. His rose inspired 
me to add the respective TI-92/V 200 treatments 
and a slider bar activity for Derive. 

Johann Wiesenbauer provides a special tool for 
DNL#60. His Sodoku-Solver is once more an ex- 
cellent example for Derive programming. 

We have an extended User Forum. A few challeng- 
ing requests asked for being answered. 

Let me finally wish you all a great 2006 - and see 
you in Dresden!! 



Don't forget DES-TIME, the DERIVE & 
CAS-TI & ACDCA-Conference is ap- 
proaching. 15 February is deadline for 
submitting papers. 

www.des-time-2006.de 


Download all DNL-DERIVE- and Tl-flles from 

http : //www . austromath . ac . at /dug/ 
http : //www . derive-europe . com/ support . asp? dug 




EDITORIAL 


D-N-L#60 


P 2 


I 


The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE & CAS- 77 User Group. 
It is published at least four times a year 
with a contents of 44 pages minimum. The 
goals of the DNL are to enable the ex- 
change of experiences made with DERIVE 
and the 77-89 /92/ Titanium! Voyage 200 as 
well as to create a group to discuss the 
possibilities of new methodical and didac- 
tical manners in teaching mathematics. 

As many of the DERIVE Users are also 
using the CAS-TIs the DNL tries to com- 
bine the applications of these modem tech- 
nologies. 


Editor: Mag. Josef Bohm 
A-3042 Wiirmla 
D'Lust 1 
Austria 

Phone/FAX: 43-(0)2275/8207 
e-mail: nojo.boehm@pgv.at 


Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English 
articles will be warmly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the 
author gives his consent for reprinting it in 
the DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE & CAS-77 Newsletter will be. 


Next issue: March 2006 

Deadline 28 February 2006 


Preview: Contributions waiting to be published 

Two Stage Least Squares, M. R. Phillips, USA 

Some simulations of Random Experiments, J. Bohm, AUT & L. Kopp, GER 

Wonderful World of Pedal Curves, J. Bohm 

Another Task for End Examination, J. Lechner, AUT 

Tools for 3D-Problems, P. Luke-Rosendahl, GER 

ANOVA with DERIVE & Tl, M. R. Phillips, USA 

Hill-Encription, J. Bohm 

Farey Sequences on the Tl, M. Lesmes-Acosta, COL 
Simulating a Graphing Calculator in DERIVE, J. Bohm, AUT 
Henon & Co, J. Bohm 

Challenges from Fermat, Bj. Felsager, DEN 
Are all Bodies falling equally fast, J. Lechner, AUT 
Modelling Traffic Density, Th. Himmelbauer, AUT 

Do you know this? Cabri & CAS on PC and Handheld, W. Wegscheider, AUT 

An Interesting Problem with a Triangle, Steiner Point, P. Luke-Rosendahl, GER 

Mathematics and Design, Hubert Weller, GER 

Diophantine Polynomials, Duncan E. McDougall, Canada 

Financial Mathematics III, M. R. Phillips 

Contour Plots and Implicit Plots, Peter Schofield, UK 

and Setif, FRA; Vermeylen, BEL; Leinbach, USA; Koller, AUT; Baumann, GER; 
Keunecke, GER and others 

Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wurmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 
Herstellung: Selbstverlag 








D-N-L#60 


Derive- and CAS-TI-User Forum / 3D-Confidence 


P 3 


Peter Wehling, Dierhagen 


P . Wehling6baf z . de 


Sehr geehrter Herr Bohm, 

Ihre Kontaktadresse habe ich von der Internetseite der DERIVE User Group. 

Seit zwei Wochen bin ich Nutzer von DERIVE 6, und ich flnde es - zumal als Nichtmathematiker - 
hilfreich, fesselnd und inspirierend. Meine Mathematikkenntnisse sind indes wohl eher als durch- 
schnittlich zu bezeichnen; ihr Erwerb liegt bereits gut 2 Jahrzehnte zuriick, und ich benutze sie als 
Nichtmathematiker - ich bin in der Zuchtungsforschung an Pflanzen tatig - im wesentlichen auch nicht 
regelmaBig und systematisch, sondern bedarfsbedingt. Als DUG-Mitglied ware ich daher kaum von 
Wert fur andere. 

Deshalb wage ich es, Sie auf diesem Wege direkt anzusprechen und Sie um Rat zu einem mathemati- 
schen Problem zu fragen. 

Es geht, allgemein gesprochen, um die graphische Darstellung von dreidimensionalen Variablenberei- 
chen, fur welche a <= f(x,y,z)) <= b gilt. 

Im konkreten Fall ist f(x,y,z) eine Funktion, welche LOD-Werte fur gemeinsame Maximum- 
Likelihood-Schatzwerte von x,y,z angibt, und zu bestimmen ist ein Konfidenzintervall der 
x,y,z-Schatzwerte, in der Form LODmax-1 <= f(x,y,z) <= LODmax. x,y und z sind ihrerseits auf defi- 
nierte Wertebereiche beschrankt. 

Die Frage ist also, wie die fur obige Beziehung giiltigen x/y/z-Variablenkombinationen gefunden und 
geplottet werden konnen, sei es als Flachenkorper oder als Wolke diskreter x,y,z-Punkte mit vordefi- 
nierten Schrittweiten. 

Ein Arbeitsblatt mit einigen Details ist angefiigt; fur Riickfragen stehe ich jederzeit gern zur Verfu- 
gung. 

Ich wiirde mich freuen und ware dankbar, wenn Sie oder ein Kollege mir einen Tipp geben konnten, 
der mir weiterhilft. (Ubrigens, fiber Programmierkenntnisse verfiige ich noch nicht. Ihr Buch "Pro- 
grammieren mit DERIVE” ist bereits geordert, wird aber wg. Urlaubszeit erst in zwei Wochen ausge- 
liefert.) 

Mit freundlichen GriiBen 


Peter Wehling 
D- 18347 Dierhagen 


Mr Wehling is a DERIVE "freshman" and wants to use DERIVE for presentating 
3D-confidence regions. Mr Wehling is not a mathematician but he uses math-tools for his 
work as a researcher in the field of plant breeding. So his request refers to results from his 
research work with genes. He describes the source of his problem in the following DERIVE 
file. This is too specific for me, so I didn't even try to translate the background into English. 

It is a statistic problem dealing with maximum likelihood estimation for three parameters. 

In a few words: The problem is to find the body which is given by all points (c,d,s) which give 
function values z(c,d,s) between 4.37 and 5.37. 0 < c, d < 0.5 and 0 < s < 1 . 
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The DERIVE-file 

Hi ntergrund : Es handelt sich um ein Problem aus der Genetik: Vier Genotypenkl assen fur 
zwei segregi erende Gene A und B mit gleichen erwarteten Haufigkeiten (1/4) und mit beo- 
bachteten Haufigkeiten M, N, 0, P. Die Genotypenkl assen sind Funktionen von drei Vari- 
ablen, die zu schatzen sind: c und d fur Rekombi nati on zwischen den Genen A-Q und Q-B 
im Interval 1 A-Q-B; s fur den sel ekti erenden Einfluss von Q auf die Haufigkeiten der 
A/B-Genotypenkl assen . c und d konnen jeweils zwischen 0 und 0.5, s kann zwischen 0 und 
1 variieren. 

Die Funktion fur den LOD-Wert zur Maximum-Li kel i hood-Schatzung von c, d und s bei beo- 
bachteten K1 assenhaufi gkei ten M, N, 0, P ist: 


#1: z(M. N, 0, P, c, d, s) := M-LOG 


4«(c*d + (I - c)»(l - d)*(l - s)) 


2 - s 


10 


' 4*(c*(l - d) + (1 - c)«d«(l - s)) 


N-LOG 


10 


2 - s 


0« LOG 


4 • C C 1 - c)*d + o(l - d)-(l - s)) 


2 - s 


10 


4 • C C 1 - c)«(l - d) + c*d*(l - s)) ^ 

P-L0G| . 10 

2 - s 


Mit folgendem Ausdruck habe ich den maximalen LOD-Wert ermittelt. Fur das u.g. Beispiel 
mit gegebenen M, N, 0, P ist z-max ungefahr 5. 37. Der dazu gehorige gemeinsame ML- 
Schatz-wert fur die c/d/s-Kombi nati on ist 0.33/0.07/0.84: 

#2: MAX (MAX (MAX (VECTOR (VECTOR (VECTOR (z( 5 , 13, 4, 26, c, d, s), c, 0.01, 0.5, 0.01), d, 

0.01, 0.5, 0.01), s, 0, 1, 0.01)))) 

Nun ist es wunschenswert , fur diesen gemeinsamen c/d/s-Schatzwert ein Konfi denzi nter- 
vall anzugeben. Fur meine Zwecke ausreichend ware z.B. das "LOD-l-Interval 1 " , welches 
alle c/d/s-Kombi nati onen mit LOD-Werten zwischen 4.37 und 5.37 enthalt. 

Man konnte sich zunachst damit behelfen, zwei dimensi onal e Konfi denzi nterval 1 e (z.B. fur 
c/d) uber den (schri ttwei se abgetasteten) Vari ati onsberei ch der verbliebenen dritten 
Variablen zu bestimmen und graphisch mit 2D darzustellen al s Menge der Tei 1 funkti onen , 
welche die betr. BOOLE'sche Verknupfung erfullen. Das habe ich versucht, und es klappt 
gut (wenn auch ein bisschen Programmi erkenntni s das schri ttwei se Abtasten der dritten 
Variablen wesentlich eleganter gestalten konnte): 

#3: (z(5 , 13, 4. 26, c, d, 1) > 4.37 A0<c<0.5A0<d< 0.5) v (z(5, 13, 4, 26, 

c, d, 0.99) > 4.37 a 0 < c < 0.5 a 0 < d < 0.5) v u*s*w 

Bei dem Versuch, auch das dreidimensionale c/d/s-Konfi denzi nterval 1 graphisch darzu- 
stellen - eine uber Schri ttwei ten regulierbare Punktewolke im 3D-Koordi natensystem wur- 
de ja schon genugen -, bin ich dann aber gescheitert. Zwar konnte ich die c/d/s- 
Kombi nati onen , fur welche die betr. BOOLE'sche Verknupfung wahr (oder falsch) ist, i n- 
di rekt mit Hilfe von MAP_LIST oder TABLE darstellen (wenn auch nicht in einem Rutsch 
fur die u.g. Schri ttwei ten und Werteberei che , weil dann mein PC wg. mangel nden Spei- 
chers streikt): 

#4: MAP_LIST(MAP_LIST(MAP_LIST(4. 37 < z(5, 13, 4, 26, c, d, s) < 5.37, c, 0.001, 0.5, 

0.001), d, 0.001, 0.5, 0.001), s, 0, 1, 0.01) 

Nur, wie konnen die betr. c/d/s-Kombi nati onen , fur welche 4.37 < z < 5,37 "true" 
ist, in Punktevektoren fur eine 3D-Graphik ausgegeben werden??? 
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DNL: You can imagine that this problem was a real challenge for me and it took me one e- 
vening to find a first - more or less - satisfying solution: 

Lieber Herr Wehling, Dear Mr Wehling, 

I tried to apply some tricks: 

I wrote a program pkte ( lower bound, upperbound) containing loops for values c, d and s. If a 
triple of co-ordinates fulfills the condition for "function value" z the respective point will be added to a 
list. 

Taking your parameters (small increments) calculation needs a long time and can lead to memory 
problems. It might be possible to extend the DERIVE memory, but I believe that the result would not 
be much better than mine with larger increments. 

The increments are given within the program using in the assigments i:+..., j:+... and k:+.... The left 
graph was created using smaller increments the cloud of points is more compact. 

pkte(u, o, i, j, k, liste) := 

Prog 

liste := [] 

k := 0 
Loop 
j := 0 
Loop 
i := 0 
Loop 

#6: If u < z(5 , 13, 4, 26, i, 

liste := APPEND([[[i , j 
WRITE([ i , j, k]) 
i :+ 0.02 
If i > 0.5 exit 
j :+ 0.02 
If j > 0.5 exit 
k :+ 0.02 
If k > 1 exit 
RETURN liste 


j, k) < o 

, k]]], liste) 


#7: pkte(4. 37 , 5.37) 



The WRiTE-command shows progress of the calculation in the status line (bottom left). Point Size 
should be set to small (issue insert Plot in the 2D-Plot Window). 
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I was not really satisfied with my result. The points appeared too large, I wanted to distin- 
guish the single points in the cloud. Peter Schofield gave an advice: don't plot a single point, 
but plot a segment with length zero. 

(By the way I remembered Euklid's classic definition of a point ) 


I had only to exchange one line in the program from above to receive the points as dots: 


Loop 

If u < z(5, 13, 4, 26, i, j, k) < o 

lists := APPENDOiste, [[[i, j, k; i, j, k]]]) 
WRITE ([i , j, k]) 


pkte_sinaTI(4. 37 , 5.37) 



Two days later I sent another mail to Mr Wehling: 

Lieber Herr Wehling, Dear Mr Wehling, 

As announced in my last mail, I was not fully satisfied with my point diagram. I developed two other 
models which should accomplish visualization in space: 

I plot the horizontal levelplanes in 0.02 z-increments and I plot the contour curves of the body built by 
the points. 

It would not be too difficult to include all parameters (boundaries and increments for c, d and s) into 
the variable list. In the existing form of the program one has to change the respective values in the 
program which is not very user friendly and might be too difficult for the unexperienced. 

With my best regards 
Josef Bohm 

If you agree I'd like to publish the result of our inspiring communication in one of the next Derive 
Newsletters. 


See my final results on the next page: 
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We produce level planes with 0.02 incre- 
ments in z-direction: 

pkte_cont(4. 37 , 5.37) 



Add 1 as third parameter to obtain only the contour lines plotted (below left). 


Scanning with 0.01 -increments in y-direction we obtain an result of better accuracy (below right). 


pkte_cont(4. 37 , 5.37, 1) 


pkte_cont2(4. 37 , 5.37, 1) 



Lieber Herr Bohm, 

von Ihren beiden neuen Entwicklungen bin ich begeistert; sie erleichtem die raumliche Visualisierung 
des Konfidenzbereichs ganz ungemein! 


Die Eingabe der Parameter direkt in das Programm ist iiberhaupt kein Problem fur mich (mein PC ist 
gerade dabei, noch engere Schrittweiten zu errechnen) - eher schon, was genau Ihr Programm mathe- 
matisch eigentlich tut; aber zu diesem Verstandnis werde ich mit der Zeit schon noch gelangen. 


Selbstverstandlich konnen Sie das Erreichte fur den DNL nutzen. 


Noch einmal meinen herzlichen Dank fur Ihre unschatzbare Hilfe. 


Mit herzlichen GriiBen 
Peter Wehling 


Dear Mr Bohm, 

I am really enthusiastic about your two new programs. They support 3D-visualization of confidence 
regions in a high degree. Entering parameters into the program is no problem for me (my PC is just 
busy calculating smaller increments) - it is a bit more difficult to understand how the program works, 
but I am sure that I will find it out very soon. 

You can use the results for the DNL, of course. Once more many thanks for your unvalueable support. 

With best regards 
Peter Wehling 



Josef Lechner: Various Leaves of the Rose-Curve D-N-L#60 


Verschiedene Blatter der Rosenkurve 
Various Leaves of the Rose-Curve 


Josef Lechner, Viehdorf, Austria, le j os6aon . at 

Abstract 


Ebene Kurven eignen sich besonders fur experimentell-heuristische Zugange zu Problem- 
stellungen. Dies soil am folgenden Beispiel verdeutlicht werden. 

Plane curves can be used in a meaningful way for experimental-heuristic ap- 
proaches for problem solving. This shall be demonstrated by the following 
example: 

Beispiel 1 Konstruktion von Astroide und Rosette / Construction of astroid and rose 

(mittels Dynamischer-Geometrie-Software (DGS)) 

Wir wollen mit einer Gleitstrecke experimentieren. Dazu gehen wir von folgender Vorstellung aus: 
Eine Leiter AB ( = Gleitstrecke ) sei an eine Wand angelehnt und wird an ihrem unteren Auflagepunkt 
waagrecht - zwischen den beiden Extrempositionen „stehend“ und „liegend 44 - hin- und her bewegt. 
Eine solche Konstruktion lasst sich mit einem DGS leicht bewerkstelligen. AnschlieBend konnen wir 
bereits mit Experimenten beginnen: Es werden dazu drei Punkte auf der Gleitstrecke markiert, der 
Halbierungspunkt, und die beiden Lotpunkte (siehe Abb. 1). Frage: „Auf welchen Bahnen bewegen 
sich die drei Punkte? 46 . Es ist eine sehr interessante Ubung, den Schulem zuerst die gesuchten Bahnen 
skizzieren zu lassen. Die Ausfuhrung mit DGS zeigt, dass ein Viertelkreis, eine Viertelastroide und 
eineViertelrosette entstehen. 

We start experimenting with a gliding segment. Ladder AB (= gliding segment) is moving between 
poitions “standing and lying". This construction can easily be performed using a dynamic geometry 
program. We mark three points on the segment - the middle point and the two petal points (see fig. 1). 
Question: “What is the locus of the three points?". It is a very interesting exercise for the students to 
let them sketch the loci before “asking" the DGS. The program shows the results immediately: there 
are appearing a quarter of a circle, a quarter of an astroid and a quarter of a rosette (if A and B are 
moving on the positive axes). 



10 


Abb. 1 / fig 1 Experiments with the gliding segment ( performed with Cabri and Geogebra) 
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Durch entsprechende Spiegelungskonstmktionen ergeben sich schlieBlich die vollen Kurven (Abb. 2). 
Die Bewegung der Gleitstrecke kann zudem mit Hilfe eines Schiebereglers gesteuert werden. 



Abb. 2 / fig 2 The moved segment delivers astroid, circle and 4-leaved rose 


Der Einsatz von Technologien (wie DGS) verstarkt experimentelles Arbeiten im Sinne von „Was pas- 
siert, wenn ... ?“ Es begiinstigt das „Fragen-Stellen“ und „In-Frage-stellen“. („Warum ist das so?“, „Ist 
das die einzige Moglichkeit?“, „Ist das die giinstigste Losung?“, ,,Stimmt das immer?“) Technologi- 
sche Werkzeuge begiinstigen das experimentelle Arbeiten besonders deswegen, weil mit konkreten 
Objekten gearbeitet werden kann, weil Bilder Assoziationen hervorrufen und bei der Entwicklung von 
Grundvorstellungen helfen. 

Experimentelle Mathematik ist Mathematik im Vorfeld einer systematischen Fachbetrachtung. Man 
sucht eine Vermutung durch systematisches Probieren. Es kann vereinfachend (nicht verfalschend) 
gearbeitet und gesprochen werden. Die Arbeit am konkreten Bespiel ist wichtiger als der Drang zur 
Verallgemeinerung, das unmittelbare Testen einer Vermutung (mit Hilfe der zur Verfiigung stehenden 
Werkzeuge) ist wichtiger als der exakte Beweis, Intuitionen und praformales Vorgehen wichtiger als 
mathematische Strenge. 

Es wird seitens der Didaktik oft beklagt, dass die experimentelle Seite im Unterricht oft viel zu kurz 
kommt; dass viel zu rasch fertige Modelle, Ableitungen, Verfahren den Schiilern prasentiert werden, 
ohne dass diesen die Moglichkeit gelassen wird, selbst Entdeckungen zu machen. Technologie- 
unterstiitzter Unterricht bietet hier vielfaltige Chancen. Naturlich kann Experimentelle Mathematik 
ohne theoretische Untermauerung relativ nutzlos sein; wenn namlich die entsprechenden Begriffsbil- 
dungen, kurz die zu Grunde liegende Mathematik, nicht angeschlossen wird. Experimentelle Mathe- 
matik kann in keinem Fall die klassische ersetzen , aber sie kann z.B. zu Beweisbediirfnis anregen , 
zum Bediirfnis, die zum Problem und Vermutungsformulieren notige Sprache zu entwickeln u.a.m. 
Kurz: sie kann helfen „die Weichen richtig zu stellen“. 

Die experimentell erzeugten, geometrisch vorliegenden Kurven mathematisch-formal zu fassen soil im 
nachsten Beispiel versucht werden. Dabei soil auf die vielfaltigen Darstellungsmoglichkeiten und auf 
eine Gegemiberstellung der verschiedenen Darstellungsformen eingegangen werden. 
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Using technologies (like dynamic geometry systems) affirm experimental working in the sense of 
"What if?“. It inspires "posing questions" and "doubting" ("Why is this so?“, "Is this the only one possi- 
ble solution?", "Is this the best solution of all?", "Is this always true?") Technological tools support 
experimental working especially because of working with concrete objects, because pictures call pro- 
voke and support developing basic concepts. 

Experimental mathematics can be a initial point of systematic view of the subject. One searches for a 
conjecture by systematic experimenting. One is able to work - and to talk - in a simplified (but not in a 
falsified way). Working on the given concrete problem ist more important than the wish for generalisa- 
tion, immediate checking of a conjecture (using the available tools) is more important than the exact 
proof, intuition and preformal acting is more important than mathematical strictness. 

Didactic often complains that math education lacks the experimental view, that too quick ready made 
models, derivations, methods are presented to the students, without leaving them a chance to make 
their own discoveries. Technology supported teaching offers many chances. Of course, experimental 
mathematics can be useless without an appropriate theoretical base, if the underlying mathematical 
concepts are not following. Experimental mathematics can in no case substitute classical mathemat- 
ics, but it can eg incite the wish for a proof, incite to develop an appropriate language for communica- 
tion about the findings a.o. To make it short: it can help to "shift the right switch". 

In the following we will show various representations of the rose curve (4-leaved rose). This is not only 
a "Window Shuttle" between algebraic and graphic representation, but also a "Shuttle" between differ- 
ent algebraic representation forms (parameter form, polar form, implicit form, 3D-representation, func- 
tion form and complex form). We hope that the mathematics behind is self explanatory, so we will 
minimize additional translations. 


Beispiel 2 Die verschiedenen Darstellungsformen der Rosenkurve 

The various representation forms of the rosette (4-leaved rose) 

a) Parameterform (parameter form) 

f X - x(0 

Definition : Wird eine Kurve k in der Form k .... x(t) = < mit t e \t x ,t 2 \ angeschrieben, 

[y=y(t) 

wobei x(t) und y{t) (stetige) Funktionen des Parameters t sind, so spricht man von der Parameter- 
darstellung der (stetigen) Kurve k. Der Parameter t reprasentiert dabei zumeist die (stetig) ver- 
flieBende Zeit oder den Winkel, der uberstrichen wird. 



Abb. 3 / fig 3 Sketch for deriving a parameter form 
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Die Parameterdarstellung ergibt sich aus folgenden geometrischen Uberlegungen (Abb. 3): 

PQ 

Fiir die Streckenlange BO ergibt sich aus cos q> = : BO = c cos cp und weiter fur die Lange r aus 

c 

r . . c 

sin m = und weiter r = BO sin<z > = c cos® sin® = — sin(2®). 

BO 2 

c 

Damit lasst sich fiir die x-Koordinate die folgende Beziehung angeben: x = r cos q> = — sin(2<^) cos < 2 ? . 

.. c 

Mit analogen Uberlegungen erhalt man fiir die y-Koordinate: y = r sin cp = — sin(2<^) sin cp 

Damit konnen wir nun in Parameterdarstellung anschreiben: 

Q 

x{cp) = — sin(2<£>) cos cp 

Rosette: x(cp)-< ^ mit <£>e[0,2;r] 

y{(p) = -sm(2(p) sin^ 


Wenn wir fiir c = 2 wahlen, so ergibt sich: 


\ Iqa 


*1: [SIN(2 $) COS(f), SIN(2 $) SIN <<|>> ] 


-1.5 -1 -Ck5 


0x5 1 1.5 


How to do it on the Tl handheld: 


Abb. 4 / fig 4 Rose in parameter form (Derive) 


jldt; 

-'xt 1 =s i n( 2 ■ t) ■ cost! t) 
-'y t 1 =s i n( 2 ■ t) ■ s i n( t) 
xt2^ 
yi2= 
xi3= 
yi3= 
xi4= 
yt4= 
xt5= 
yt-5= 


tep=nx5G 
in= -2.975 


cl = l. 
in= "1 . 27 
ax= 1.275 
cl = l. 


I MAIN ROD flUTD F hF; 


Abb. 5 / fig 5 Rose in parameter form (TI) 


b) Polardarstellung (polar form) 

Definition: Wird eine Kurve k in der Form 

r = r(tp) mit cp&[(p v (p 2 \ 

angeschrieben, so spricht man von der Polardarstellung der Kurve k. Polardarstellungen eignen 
sich besonders fiir Kurven, die symmetrisch um einen Punkt (meist den Ursprung) liegen. 

Die Polardarstellung ergibt sich sehr einfach unter Verwendung der oben hergeleiteten Parameterdar- 
stellung: 

c 2 c 2 c 2 

r 2 = x 2 + y 2 = — sin 2 (lep) cos 2 cp H sin 2 ( 2cp ) sin 2 cp = — sin 2 ( 2cp ) . 
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Rosette: 

Q 

r{(p) = — sin(2<£>) mit cp e [0, 2n\ 


2 


Wenn wir wieder c = 2 wahlen, ergibt sich die bekannte Rosette, diesmal in Polardarstellung: 


/2: i’ = SIN (2 f) 



c) implizite Darstellung (implicit form) 

Definition: Wird eine Kurve k in der Form 


f(x,y) = 0 


angeschrieben, so spricht man von der impliziten Darstellung der Kurve k. Treten in dieser Dar- 
stellung nur Potenzen und Produkte von x und y auf, so spricht man von einer algebraischen Kur- 

ve: Y^Y J a km -x k -y m = 0. 

k = 0 m = 0 

Die implizite Darstellung lasst sich problemlos aus der Polardarstellung gewinnen: 

Unter Bemcksichtigung der beiden Transformationsbeziehungen x = r cos cp und y = rsmcp erhalten 
wir 

2 2 

2 2 2 *2 2 ^ y i ’ t 2 2 2 2 

r =c cos ^sin (p = c — — — und welters ( r ) = c x y . 

r r v 7 

Da aber r 2 = x 2 + y 2 ergibt sich 

Rosette: (x 2 + y 2 f - c 2 x 2 y 2 = 0 

2 2 3 2 2 

/3: (x + y ) = 4«x «y 



Abb. 7 / fig 7 Rose in implicit form 
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On the Tl produce a contour plot of level 0 (see below). 



J'LDT:' 

■^zl=(x 2 + y : 
z2= 
z3=H 
z4= 
z5= 
z6= 
z7= 
zS= 
z9= 

2) 3 _ 4 . x 2. y 2 

s3<x,y>= 


MAIN 

RAD AUTD 5D 


fVJSIzoonl 1 

eyey=jB^DH 


eye+=EL 
eyeT=0. 
xni i n= - 1 . 
xnax= 1 . 
xgrid=2G. 
yr-i i n= - 1 . 
ymax=l . 
ygrid=2G. 
zni i n= - 1 . 
zmax=l . 
nconiour=l , 


MAIN 

RAD AUTD 3D 



Or do immediately an implicit plot. 




Abb. 8 / fig 8 Rose in implicit form on the TI 


d) Darstellung als Funktion in zwei Variablen (function of two variables) 


Die implizite Darstellung einer Kurve k lasst sich interpretieren als Schnitt einer Funktion z(x,y ) 
in 2 Variablen mit der Ebene z = 0. 

(x 2 +y 2 f-c 2 x 2 y 2 =0 

z (x,y ) 


Die Kurve erscheint in dieser Betrachtungsweise als Hohenlinie (die Kurve bekommt damit sozusagen 
ein „dreidimensionales Leben 44 ). Mit CAS-Unterstiitzung lasst sich dies gut visualisieren. Die Schuler 
konnen diese Darstellung auch leicht drehen und aus verschiedenen Blickwinkeln beobachten. 



Abb. 9 / fig 9 Function of two variables - 3D plot 
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The right sequence of screeshots show an attempt to 
visualize the surface together with the intersecting plane 
z = 0. 



JLOTl- 

= l= C 




x 2 + u 2) j _4. x 2. u 2 
,y),zl(x,y)>0 


^z3=naxUx 

z4 : 
z5 ; 


((x 2 + y 2 ) 


■ 4 ■ x J 





xnax=l . 
xgrid=35. 
mm i n= " 1 . 
ynax=l . 
ygrid=35. 
znin= 01 
znax=.01 
ncont-ourm 


I'FiTTiYTF^T - F5 Y FH YT?^YTI^tf7 
f— |ZQQn|Trace|Regraph|Maih|Dratj|T 



MAIN RAD AUTD 3D 


Abb. 10 / fig 10 Function of two variables in Derive and on the TI 


DERIVE verfugt iiber keine implementierte Funktion zur Darstellung von Konturkurven (Hohen- 
schichtenlinien). Auch in den USER-Dateien kann man nichts derartiges finden. Ein Artikel im na- 
achsten DNL wird helfen, diese Liicke zu schlieBen. Abbildung 1 1 zeigt einige Konturkurven mit der 
extra betonten 0-Kurve (links) und die Flache mit der O-Niveau-Kurve (rechts). 

Derive does not offer an implemented Contour Plot function. There is no such function among the 
USER files. Refer to a contribution in the next DNL and learn how to fill this gap. Figure 11 shows 
some level curves with the 0-level curve (left) and the surface with the 0-level curve (right). 




Abb. 11 / fig 11 Contour plot with Derive 
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I failed trying to perform the same procedure on 
the TIs because the handheld device is unable to 
solve cubic equation. What we can do is to trans- 
fer the DERIVE solutions to the Tl and then plot 
the rosette as composition of functions. 

See below Figure 13 


I send expressions #7 and #8 from above (using solutions instead of solve) to the V200. In the Tl- 
Text editor I accomplish the script storing the elements of the list of solutions as functions yl(x) 
through y6(x). Then I execute the script. Inspecting the table I find that two branches are containing 
only imaginary points, so they can be omitted for the plotting procedure. (By the way, try plotting ex- 
pression #8 from above with Option "Plot real and imaginary parts" activated - Derive, of- 
course!) The screen shots below illustrate how it was done: 


i'fIT^t FF? YTF^t fh t y fl t Fe^ t "i 
IflIgebralCalcloiherlPrgnlolClean Up| [ 


■ £olue((x 2 + y 2 ) 3 = 4 ■ x 2 ■ y 2 , y) 

y^(y^ + 3y 2 x 2 + x 2 (3x 2 -4))= - x 

■ solue((u + = 4 ■ u ■ w , w) 

w ■ ( w 2 + 3 ■ u ■ u + u ■ ( 3 ■ u - 4)) = ■ u ^1 


solue<<u+w> A 3=4u*w, 


IHqrninj: Mort solutions roqy ■Lxist 


Abb. 13 / fig 13 Try with the Tl 




-YT^rY ™ T - ft Y 1 

J Uiew Execute Find... 

... X A T _ 4 * , ! A V , , 


U : t2eros ( i ; x x 2+'=r2 ; A 3- ; 4*x A 2*y A 2 ; i , y ; 

> -i < 3 !>*x > V3 , -3 1 * < 3x4 < x* < 4*£ i n < s i n* < 
3*f(3)+x/4)/3)-^3)^) >x3,3 A <3x4>*r< -x 
*<2*J‘<3)*sign<x)*cos<sirvi<3*J'<3)*x^4)^ 
3 > +2*s i n < s i n* < 3*J < 3 !>*xx4 V3 > +T < 3 !>*x > !>x 
3, -3-- C 3--4 > =+=-4" C -x*C2*T<3:>*£ign<x:>*c0£(:£i 
n-i < 3*T < 3 ) *xx4 ) x3 ) +2*£ i n < £ i n* < 3*J < 3 ) *xx 
4 > x3 > +T < 3 !>*x >> x3 , 3' * < 3x4 < x* < 2*T < 3 !>*£ 

ignCx>+c0£(:£in-i(:3+-J'(:3>+xx4>x3>-2+£in(:£ 
in*<3*r<3>*xx4>x3W<3>*x>>x3, -3 A <3x4> 
< x* ( 2*J (3>*£ign(x> *co£ < £ i it 1 < 3*T ( 3 > *x 




'y5= 


-2-J3 ■ £ign(x) ■ co£^ ^"3” ► 


'y6= 

_aZz_ 


-2 ■ J3 ■ £ign(x) ■ co£ 




v6<x>=...<3>ttx/4>/3>+J<3>ttx>>/3 


[t ^^jcomandju i eujE xecute |F i nd. 


x) yx3J+£ol 
CS£Ol [1 ] 

C: an£< 1 )+yl(x) 
C: £ol [2] 
Csan£<l)+y2<x) 
C: £ol [3] 

C: an£< 1 >-*y3(x> 
C: £ol [4] 
Csan£<l)+y4<x) 
C: £ol [5] 

C: an£< 1 >-*y5(x> 
C:^ol [6] 


E5L 



X 

y i y2 

y3 

y4 

y5 


5GGG 

. 1399S 14GG 


-1.16... 

. 76653 

4GGG 

. GS555 GS5G 

1 . G15... 

-1.G1... 

. 73696 

3GGG 

. G4GG4 -.0466 

. S54S... 

-■S54... 

.67718 

2GGG 

. G2G31 G2G3 

. 6772... 

-.677... 

. 5816G 

1GGG 

. GG5G2 GG5G 

. 4634... 

-.463... 

. 42988 

G.GGGG 

G.GGGG G.GGGG 

G.GGGG 

G.GGGG 

Q.QQQQ 

. 1GGGG 

. GG5G2 -. GG5G 

. 4634... 

-.463... 

. 42988 

. 2GGGG 

. G2G31 -. G2G3 

. 6772... 

-.677... 

. 5316Q 

v3<x>=1.165*i 

MAIN RAD AUTD FUNC 


PiSblffir.IcK’cWSi-lPrilol 

Clean UrJ I 

-3 3 ' 4 - 


_ - r= . . . 

t. 4 

X 

jl -J ■— 1 — ■ 1 91^ Ik XJ cos 

3 > 




3 

- 3 3X4. 


ji. ■ -J ■ si9ri( x) ■ cos 

C ■ J 3-J3 

sin H— 4 . 

X 


L 3 

3 

Done 


.<3>*x/'4>/'3>+J<3>*x>>/3-»-j6<x>I 


func ed/so 




> 




Abb. 14 / fig 14 InterConnectivity between Derive and the Voyage 200 


f) Darstellung in der komplexen Zahlenebene (complex form) 

Die Rosenkurve besitzt - kaum iiberraschend - auch eine Darstellung im Komplexen. Diese er- 
halten wir, wenn wir von der Parameterdarstellung ausgehen und die beiden Koordinatenfunktio- 
nen als Real- und Imaginarteil interpretieren. 


x(t ) = sin(2f) cos(0 
y(t) = sin(2f) sin(f) 


z(t) = sin(2^) (cos(0 + i sin(f)) 


Unter Verwendung des EULERSCHEN Satzes erhalten wir schlieBlich: 


Rosette: z(t) = sin(2f) e lt 


Mit einem kleinen Zeichenmodul „cplot ( z_) das wir uns definieren, erhalten wir die Rosette in 
der GAUSZschen Zahlenebene. 

#9: cp1ot(z_> :z [RE(z_) , IFI(z_)] 

#10: z_rosette(^> :z SIN (2 COS(^) + i SIN (2 f) SIN(f) 

#11: c plot (z_rosette (^)) 









Abb. 16 / fig 16 Complex form on the T1 


Weitere Aufgaben konnten gestellt werden, wie zB: 

• Lasst sich die Astroide auf ahnliche Weise so erschopfend behandeln? 

• Welches ist die Ortslinie eines beliebigen - festen - Punktes auf dem gleitenden Segment? (Da- 
mit gelangen wir zur Ellipsenbewegung und zur punktweisen Konstruktion einer Ellipse = Pa- 
pierstreifenkonstruktion). 

• Erzeugung der Rosenkurve und der Astroide als Radlinie. 

Follow up tasks could be: 

• Can you treat the Astroid in a similar way? 

• What is the locus of any fixed point of the gliding segment? (We find the ellipse motion and a a 
way to sketch an ellipse point by point.) 

• Create the rose curve and the astroid as a trochoid. 











I was vey much inspired by Josef's contribution, so that I did not only add the Tl-solutions - 
with his agreement, of course, but had additional ideas to open another “Presentation 
WINDOW”. 


Mit Hilfe der Schieberegler lasst sich mit DERIVE eine Liicke zwischen dynamischer Geometrie und 
Computeralgebra zumindest teilweise schlieBen. Die Konstruktion unter Verwendung von Parametern 
lasst einerseits die Animation der Grafik zu, zwingt aber andererseits sofort zu einer notwendigen 
Verallgemeinerung der „konkreten Objekte”. Die angestellten Uberlegungen werden sofort verifiziert 
oder falsifiziert. Dies eroffnet eine ganz neue Dimension im Mathematikunterricht. 

Presentation of the "gliding ladder" and the interesting three points using two slider bars helps to close 
the gap between dynamic geometry and computeralgebra. I insert two slider bars, one for the length r 
of the ladder (with 0 < r < 6) and the other one to move the base point A(f,0) of the ladder (with -6 <t< 
6). We calculate and plot simultaneously in two open windows (Algebra and 2D-Plot). Parallel geoem- 
tric and generalized algebraic representation opens new dimensions in math teaching. 



Midpoint and the intersection points of the perpendicular lines with the segment are calculated - and 
plotted. 


The co-ordinates of midpoint m represent the 
parameter form of the locus of m. 

For plotting the locus we have to change the 
variable for the parameter, because variable t 
is now preserved for the slider bar. 


1 

t&: ii := a + ■ (b - a) 

2 






2 



(Are the students able to recognize the curve 
from its parameter form?) 


# 10 : 


2 


2 
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#11: SOLVE (ladder = perpl, [p, t| ] ) 


/ 12 : 


2 2 11 
i- - t t -J(r - t ) 

p = * q = 


#16: SOLVE (ladder = perp2, [p, t|]) 

1 

/17: 


1 1 

t t 7(r - t ) 

P = * q = 

1 1 


#13: SUBST 


#14: 


r 2 2 ” 


r 2 ' 

r — t 


t 

ladder, p, 

/IS: SUBST 

ladder, p, 

2 


2 

l i' j 


^ r j 


3 2 23/2 

t (r - t ) 


/IS: 


3 2 23/2 

t_ (r - t_ ) 


# 10 : 


/ 20 : 


3 2 2 2 

t t -7(r - t ) 


t - 


3 2 2 2 

t_ t_ 7(r - t_ ) 


The plot shows the three loci - but not as a static result but as a "living" object of our own creation. 
By the way we found another parameter form for the three curves (without trig functions - only using 
analytic geometry!). 




Abb. 18 / fig 18 All loci and slider bars for r and t 
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#21 creates the whole curves in one step: 


# 21 : 




Abb. 19 / fig 19 The parameter forms and the respective plots 


Eine nette Aufgabe ist es auch, aus dieser Parameterdarstellung die implizite Form der Rosenkurve zu 
gewinnen. Zuerst wird uberpriift, ob die Parameterdarstellung in die vorliegende implizite Form passt. 
Dann ist es wichtig zu erkennen, dass auch in der Komponenten fur x der Term r 2 -t 2 enthalten ist. 


Einige kleine Schritte werden im Kopf- oder am Papier - durchgefiihrt, die restliche Prozedur erfor- 
dert auch eine gehorige Kompetenz im Umgang mit dem Programm, aber auch im Erkennen der not- 
wendigen Strategic, nur Verrechnen und Schreibfehler sind Geschichte. 


It is a nice exercise to derive the implicit form of the curve from this parameter form. At first we check if 
the parameter form "fits" into the given implicit expression. Then it is important to recognize that ex- 
pression r 2 - f is hidden in the first component of the parameter form. 

Some minor steps are performed mentally, the remaining procedure requires a good deal of compe- 
tence using the tool, but also deciding the right strategy. But Calculating and writing errors are "histo- 
ry". 


SUBST 


2 2 3 

Of + Y ) 


2 2 2 

r x y , [x, y], 


3 2 2 2 

t t - t ) 


t - 


2 3 

r t - t 


/40: 


true 


/4i : 


2 2 2 2 2 
t (r - t ) t f (r - t ) 

x = , y = 


/ 42 : 


2 2 2 1 
t 4(v ~ t ) 


4 2 2 

t (r - t ) 
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#43: 


t w 2 t w 

sf = , y = 


#44 : SUBST 


2 t w 


w, 


i' 


3 ^1 

2 t x 


#45 : SUBST 


2 2 
t Ci’ - t ) 
x = , t, 


"22" 

y i’ 

1/3' 


2/3 2/3 2 

SICNCx ) |r y| y 

v — 


J 


A — 

1/3 x 

X j 


#46: 


2/3 2/3 2 

SICNCx ) 1 1 ’ y| y 


x = 


1/3 


2 2/3 2/3 'l 

y SICN(x ) 1 1- y | 


x + 


1/3 


#47: 


2 2/3 2/3 

y SICNCx ) 1 1 ’ y| 


x + 


1/3 


2 2 3 2 2 ^ 

Cx + y ) r y 


#48: 


2 2 3 2 2 

Cx + y ) r y 


#49: 


We did it! 

2 2 3 2 2 2 

(x + y ) = r x y 


It might be easier for students to work with r = 2. 

See finally the Cabri-version performed on the Tl: 





. 

y 7 

\ 1 
V i 
_A i 



1 

1 

1 

MAIN 3 DEG EXACT FUNC 


I- 1 * 1- ■ ✓ l?i2 1 

^ ^ 1 
■v 1 

jj J 


X ! 

■■-j ^ ^ 

yT/\x. h 

.-■■■ , T^_ 

f-'. 

jr 

1 

1 

1 

■v V 

MAIN 3 DEG EXACT FUNC 


Abb. 20 / fig 20 The Tl-Cabri version 
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Numerics versus Symbolics 

B Kutzler (Linz, Austria) 

This lecture is a meditation about two concepts which, in the context of computer 
algebra, sometimes appear as opposing each other. 


1. Etymology 

We start with the etymology of the two words. The word “numerics” comes from the Latin word “nu- 
merus” which means “part”, “number”, where “number” itself is derived from “part” as the result of a 
counting process. 

The word “symbolics” comes from the Greek word “symbolon” which is composed of the two words 
“sym” (meaning “together”) and “ballein” (meaning “to throw”). Therefore, “symbolics” means “to 
put together”. “Putting together” can be for two reasons: It can be for constructing something, i.e. for 
creating a whole from parts. And it can be for putting things next to each other so that they can be 
compared 1 . For the Greek a “symbolon” was anything that would be comparable to the real thing 
whose place it took. 

It is important to note that, strictly speaking, “5” also is a symbol for, say, the number of fingers of a 
hand, and “3.2” is a symbolic representation of the number obtained by dividing 32 by 10. But we 
don’t use “symbolics” in this narrow sense of the word here, because then all mathematics would have 
to be called symbolic mathematics. 2 

The German mathematician C F Gauss said: “Mathematics is concerned only with the enumeration 
and comparison of relations. ” With the above in mind, this means that mathematics is concerned only 
with numerics (“enumerate”) and symbolics (“compare”). 


2. First thoughts 

In a computer algebra system (we use Derive 6) enter yp2A , simplify, then approximate. 


# 1 : 

42A 


#2: 


2 -Vs 

# 3 : 


4 . 89 89 794 85 


#1 and #2 are two different symbolic representations of this number. #3 is a numeric (decimal) repre- 
sentation of a ten-digit approximation of the same number 3 . Pragmatically, we can say that numeric 
mathematics is the mathematics on numbers (typically in decimal notation) such as #3. Symbolic 
mathematics is mathematics on everything else. 


1 The word “compare” is composed of the two Latin words “com” = ”cum” (meaning “together”) and “par” 
(meaning “equal”). 

2 In fact, in the narrow sense of the word “symbolic” one could consider mathematics the science of symbols. 

3 Compare our respective comment from section 1 . 







D 


-N-L#60 


Bernhard Kutzler: Numerics versus Symbolics 


p23 


With a calculation such as the above on the screen of a symbolic calculator, Bert Waits once asked: 
“How do you recognize a mathematician? ” and suggested the following answer: “A mathematician 
considers #2 a beautiful result. ” So we mathematicians like symbolics more than numerics ... al- 
though we know that numerics has its virtues too, and sometimes we cannot do without numerics. 
More about this later. For now we quote C F Gauss again, who has said that “a poor mathematical 
education often is demonstrated by a highly developed skill of mental arithmetic. ” (In the context of 
this paper this could be rephrased as: “Poor symbolics often is demonstrated by good numerics. ”) 


As we said, expression #3 represents only an approximation of n/ 24 . The precise decimal representa- 
tion of this number has infinitely many digits and, therefore, cannot be written in a finite amount of 
time or in a universe with only a finite amount of matter. Therefore, practical numeric mathematics 
necessarily is an approximative mathematics. 


3. Different kinds of mathematics 

The integers I, the rational numbers Q, and the real numbers R are three important number sets in 
(school) mathematics. They possess useful properties. An important property is “closure”. It guaran- 
tees that one remains inside the domain when performing an arithmetic operations. I is closed w.r.t. 
addition, subtraction, and multiplication. Q and R are closed w.r.t. addition, subtraction, multiplica- 
tion, and division (except for 0). The “deficiency” of the integers is that division can take us outside 
the domain, for example when dividing 3 by 4. 

In Q (and R) we can divide 3 by 4. The result, % or 0.75, is an element of Q (and R). If you divide 1 
by 3, the result, 1/3 or 0.33333..., also is an element of Q (and R). But the latter example causes “trou- 
ble” when it comes to a “material” representation of the number in “practical” numeric mathematics, 
for example on a computer or a calculator, where numbers are represented in decimal notation with up 
to n digits (sum of digits before and after the comma). On calculators n often is 12 or 14. 

On such a device, the result of the division of 1 by 3 cannot be represented, therefore an approxima- 
tion obtained by cutting off (infinitely many) decimal places is used. For example 1/3 will be ap- 
proximated by 0.33333333333333. 

Floating point arithmetic uses decimal numbers with up to n digits total. Fix point arithmetic uses 
decimal numbers with up to n digits following the comma. Let’s look at the somewhat simpler case of 
fix point arithmetic. We formalized this as follows: Define R(n) to be the set of all real numbers with n 
decimal places. Then: 

I R(0) d R(l) ci R(2) c= ... d R (n) d ... d Q d R 

I is closed w.r.t. addition, subtraction, and multiplication 
R(/i) is closed w.r.t. addition and subtraction 

Q and R are closed w.r.t. addition, subtraction, multiplication, and division 

Closure is an important property when it comes to the validity of identities. For example, the simple 

identity — -x = 1 is valid in Q and R, but it is not valid in any R(^). On most calculators this fact is 
x 

hidden in “obvious” cases because numbers such as 0.999999999 are rounded to 1. What really hap- 
pens can easily be visualized in Derive. The following is a graph of y = — -x — 1 in R(l). 

x 
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As you can see, the graph mostly is different from 0, hence for most values of x — • x ^ 1 . Screen im- 

x 

ages for R(2) and R(3) follow. Appendix 1 gives the code of the function used to produce these 
graphs. 



■0.05 

X 


^ : :e ' : :: : : A 

■^0.05 

•l-o.l . . . . 



Another example is the identity (x 2 - y 2 ) = (x + y) • (x - y) . In R(l) we have 


(1 . 1 2 - 0.2 2 ) = 1 .2(1) - 0.0(4) = 1 .2 

and 

(1.1 + 0.2) •(!.!- 0.2) = 1 .3 • 0.9 = 1 . 1(7) . 


#S: 

testCl.l, 0.2) 


#9: 


[1.2, 1.1] 


Also for this example, falsity of the identity can be visualized in Derive. The following picture shows 
the graph of z = (x 2 - y 2 ) - (x + y) • (x - y) in R(l). The result is different from zero for many val- 
ues ofx and y. 
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With the help of Derive’s trace function one can easily find concrete examples in any R(n). We find 
x = 1 , y — 0.5 as another example in R(l) (above picture) and x = 0.02 , y - 0.0075 in R(4) (below 
picture). 



None of the traditional computation tools such as calculators are suitable for the “classical” mathemat- 
ics in Q or R. Only computer algebra systems with their symbolic representations of rational numbers 
as quotients of two integers, with fractional powers, with re, e , etc., are appropriate. 

By the way, the traditional tool for numerics is the abacus. Slide rules, four-function-calculators, sci- 
entific calculators and traditional (numeric) computer software are but sophisticated editions of an 
abacus. Computer algebra systems are a quantum leap. They are for symbolics what the abacus is for 
numerics. 

It is perfectly fine to do mathematics in an “R(n) environment”, but students need to understand the 
consequences. They need to understand, ideally by experiencing it with appropriate examples, what 
can happen with identities such as the above. In Derive one can do approximate arithmetic with a 
specified number of digits. This is of great help when studying or demonstrating the effects in an 
“R(n) environment”. 
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We end this section by looking at the example x ■ \[x ■ (y[x + 1 + V x - 1 - 2\[x j . Enter the expression, 
substitute 1 million (1 000 000) for x , then approximate. 


#i: x.7 x, (7( x + i) + 7o - i) - 2«7 >o 

#2: 1000000 -71000000 .(7(1000000 + 1) + 7(1000000 - 1) - 2-71000000) 

#3: 0 


Looks like a clear zero. By default, Derive uses ten digits for approximations. If, instead, we approxi- 
mate #2 with an accuracy of 15 digits, the result is very much different: 


#4: -0.250012501243889 


In [Kutzler/Kokol-Voljc 2003], pages 72ff, we show how to use Derive 6 to investigate this example 
in detail and demonstrate how it happens that the two results are so different. 


4. Numerics helps Symbolics (i): Limitations of Symbolics 

Symbolics has its limitations. Some symbolic computations are impossible, 

- because an algorithm cannot be found due to theoretic limitations, 

- because an algorithm has not been found yet, 

- because an implementation of the algorithm does not exists yet, 

- because the execution of the algorithm requires too much time or space. 

In some of these cases a numeric solution is better than no solution. 

Examples are: 

(a) It is impossible to find symbolic solutions (using known functions and constants) of general poly- 
nomial equations of degree higher than four. 


5 4 3 2 

#4: x + 2 ■ x - 3 ■ x + 4 ■ x - 5 ■ x + 6 = 0 

5 4 3 2 

#5: S0LVE(x + 2-x - 3-x + 4-x - 5-x + 6 = 0, x) 

5 4 3 2 

#6: x +2-x -3-x + 4 ■ x - 5 ■ x = -6 

5 4 3 2 

#7: NSGLVE (x + 2-x - 3-x + 4-x - 5-x = -6, x) 

#8: x = -0.2389321333 - 1 . 151781933 - 1 v x = -0.2389321333 + 

1 . 151781933 ■ l v x = 0.9390814243 - 0 . 6271970073 - i v x = 
0.9390814243 + 0 . 6271970073 ■ l v x = -3.400298582 


D 


-N-L#60 


Bernhard Kutzler: Numerics versus Symbolics 


p27 


(b) For certain functions it is impossible to find closed form antiderivatives (using known functions 
and constants). 



5 





i — i 



SIN 


dx 

m 





0 



#5: 

2 . 035618287 


(c) It is impossible to find, for certain classes of expressions, an algorithm which can decide the 
equivalence of two expressions. 

(d) Constructng a Groebner basis has double exponential complexity. 


5. Numerics helps Symbolics (ii): The Creativity Spiral 

When asked how he came upon his theorems, C F Gauss answered: "... through systematic, palpable 
experimentation. ” 

According to one of the epistemologically oriented theories one can visualize the main steps of 
(mathematical) discoveries as follows: Applying known algorithms produces examples. From the ex- 
amples we observe properties, which are expressed as a conjecture. Proving the conjecture yields a 
theorem, i.e. guaranteed knowledge. The theorem 4 s algorithmically usable knowledge is implemented 
in a new algorithm. Then the algorithm is applied to new data, yielding new examples, which lead to 
new observations, ... 


Apply 

opply 



Find(Experiment) 


new examples 
examples 


algorithm 
new algorithm 



implement 


theorem 

Secure 



observe 


conjecture 


i prove 


This picture of a spiral which demonstrates the path of discovery of (mathematical) knowl- 
edge was proposed by Bruno Buchberger. A detailed description of Buchberger ’s Creativity 
Spiral and references to related models can be found in the highly recommended (German 
language) book [Heugl/Klinger/Lechner 1996]. 
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In this spiral we find three phases. During the phase of finding/experimenting one uses known 
algorithms to generate examples, then obtains conjectures through observation. During the 
phase of securing conjectures are turned into theorems through the method of proving, then 
algorithmically useful knowledge is implemented as algorithms. During the phase of applying 
one applies algorithms to new data. 

In the finding phase often we have to find patterns in sequences of numbers. A well known example 

n 

from the topic “proof by induction” is to find a closed form expression for For this particular 

i=\ 

example there exists an elegant solution which is based on an observation made by C F Gauss when he 
was still very young. Here we use a method which may be helpful also for other such problems. 
Clearly, there are many alternative approaches which will lead to the solution. 

We compute the sum for a few (consecutive) values of n , for example n = 5,6,7,8,9,10 . The results 
are stored in the second column of the below table. In looking for a pattern we perform a simple fac- 
torization of the sums into two “obvious” factors - and store the result in the third column. The pattern 
is striking, but still a little hard to describe. By doubling the first factor we make it closer to the second 
factor (fourth column). Now we see that we always have products of two consecutive numbers. In the 
fifth column the two factors are arranged in ascending order, in the sixth column we divide by two 
again to compensate for the earlier doubling. 


n 

n 

I' 

i = 1 

simple fac- 
torization 

double the 
first factor 

sort fac- 
tors 

divide by 
two 

5 

15 

3-5 

6-5 

5-6 

(5 • 6) / 2 

6 

21 

3-7 

6-7 

6-7 

(6 • 7) / 2 

7 

28 

4-7 

8-7 

7-8 

(7 - 8) / 2 

8 

36 

4-9 

8-9 

8-9 

(8 • 9) / 2 

9 

45 

5-9 

10-9 

9-10 

(9 - 1 0) / 2 

10 

55 

511 

10-11 

10-11 

(10 -11)/ 2 


Now the pattern is obvious (at least for these 6 values of n) and we can come up with the conjecture 
^ n • (n + 1) 

/ i — , which can easily be proved by induction. 

i = 1 2 

This is an inductive process which leads from the special to the general, i.e. from numeric to symbolic. 

The other two phases in the above spiral, securing and applying, are both deductive processes leading 
from the general to the specific. The applying phase is the reverse of the finding phase’s inductive 
process: You go from symbolic to numeric, for example by using the closed form expression to com- 

100 100*101 

pute the sum of the first 100 natural numbers as = = 5050 (tribute to C F Gauss!). 

i= i 2 

The securing phase also is a deductive process, i.e. it leads from the general to the specific, but this is 
on the level of mathematical logic, where the general are the (inference) rules of logics and the spe- 
cific is the mathematical theory in which the conjecture is formulated. 
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6. Symbolics helps Numerics (i): Preprocessing 


Say we need to calculate the perpendicular bisector of two points a and b . The steps of using the so 
called normal vector form are simple. 


First we compute the midpoint of a and b . This point, we call it p , is a point on the line we are 


looking for. Then we compute the vector from a to b . This vector, we call it n , is normal to the line 
we are looking for. Using n and p we can write the equation of the line using the normal vector 


form(ula) n- 


= n- p . These mostly numeric computations are easily performed with paper and 


pencil, we use Derive: 



If we need perpendicular bisectors very often, for example because we do analytic geometry, the 
above procedure becomes tedious. It is always the same steps, only the numbers (i.e. the four coordi- 
nates of the two points) are different. 

To save us from these many numeric applications of the above procedure we can apply it one time to a 



fU 


(xA 

pair of symbolic points, i.e. points a and b with symbolic coordinates: a = 

\y a j 

and b = 

Aby 


In a paper and pencil environment this “symbolic application” of the procedure is much more work 
than its meric application (and most students don’t like this kind of calculations). In a computer alge- 
bra system the extra effort of a symbolic application is done by the machine, so it is not any more dif- 
ficult for the user. Below is the computation performed in Derive. 
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bisector of a and b . Now for any new points a and b we can simply substitute the coordinates into 
this expression and obtain the resulting equation of the line. 

The one-time investment of a symbolic application saves a potentially infinite number of numeric 
applications. This method can be considered a “preprocessing”. 

7. Symbolics helps Numerics (ii): The Pentium-Bug 

On Oct 19, 1994 Dr. Thomas R Nicely, professor of mathematics at Lynchburg College discovered 
what later has become known as the Pentium Bug. After lots of testing he found that 824,633,702,441 
divided by itself gave 0.999999996274709702 - instead of 1. Later it was found that this is true for all 
numbers between 824,633,702,418 and 824,633,702,449. 

This is the well known part of the story. Here comes the unknown part 4 : Shortly after the Pentium Bug 
was found, the authors of Derive, David Stoutemyer and Albert Rich from Soft Warehouse, Inc., sent 
a complimentary copy of Derive to Prof. Nicely. After several weeks Thomas Nicely phoned up 
David. He apologized for taking so long to respond, but after he became famous for having found the 
Pentium Bug he received special protection from the police and all mail sent to him had to undergo a 
screening by a security team. Therefore he got the package with some delay. He thanked David for 


4 Personal communication with David Stoutemyer. 
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the copy of Derive and said that, in fact, he owns a copy of Derive and that Derive ’s ability to do ap- 
proximate arithmetic to a specified number of digits (our above “R(n)-arithmetic”) has helped him to 
locate the cause of the problem in the Pentium chip. 


Therefore, the bug in the numeric methods of the Pentium chip was discovered with the help of a 
symbolic tool. 


8. More from etymology 

In section 1 we said that “numerics” comes from “numerus” (= “part” (of a whole), “number”) and 
“symbolics” comes from “symbolon” (= “to put together”, “to compare”). In mathematics we use only 
the meaning “numbers” for “numerics”. “Symbolics”, on the other hand, can denote either a composi- 
tion (of parts) or something which takes the place of something else (and, hence, is comparable to it). 

Working with decimal representations of numbers is what we call “numeric”. A variable, for example 
x, which takes the place of something, we call “symbolic”. Another example of something “symbolic” 

is y = x 2 . y = x 2 stands for a whole which is composed of (infinitely many) parts, namely the pairs 

(x,x 2 ) for x taken from a certain set. When we take x from the real numbers, these pairs can be 
interpreted as points in the plane forming a parabola. 

A famos quote from the Greek philosopher Aristoteles says: “The whole is more than the sum of its 
parts ” In the above example, the points of the parabola are the parts, the parabola is the whole. The 
parabola is composed of infinitely many points, so it is the sum of its parts - but it is more than that, 
because it has properties, which none of the points has or which none of all possible subsets of the 
points has. Such properties are continuity, symmetry, or simply the fact that the points of the parabola 

are exactly those points of the plane whose coordinates satisfy the equation y = x 2 . 

Instead of “symbolic computation” we often see “algebraic computation”. How are these two words 
related? The word “algebra” comes from the title of a work written around 825 by the Arabic mathe- 
matician known as al-Khowarizmi entitled “al-jebr w’al-muqabalah”. In Arabic “al-“ is the definite 
article “the”. The first noun in the title is “jebr”, which means “reunion of broken parts”, from the verb 
“jabara” which means “to reunite, to consolidate”. 5 The second noun in the book title is from the verb 
“qabala”, with meanings that include “to place in front of, to balance, to oppose, to set equal.” 

Together these two words describe some of the manipulations so common in algebra. 


9. Pythagoras’ View of Mathematics 

Pythagoras of Samos (580 - 496 BC) was a contemporary of Confuzius (551 - 479 BC) and Prince 
Gautama, the Buddha, (560 - 480 BC). He is considered one of the wisest men of antiquity. Pythago- 
ras lived on Samos, a Greek island. The tyrant Polycrates took power in 538 BC (first with his brother, 
then alone in 532 BC). Pythagoras disagreed with his rule and left the island 532 BC for South Italy, 
where he started a school in Croton. 


5 This corresponds to the one meaning of the Greek word “symbolon” = “to put together”. 
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Pythagoras used three disciplines of teaching: 

- Nutrition for cleaning and developing the physical body 

- Music for cleaning and developing the emotional body (the soul) 

- Mathematics for cleaning and developing the mental body (the mind) 

Pythagoras had two ideals: Freedom and philosophy. For him, mathematics was the discipline to ac- 
quire both as I will explain in the sequel. Pythagoras regarded the three bodies (physis, emotio, ratio) 
to be closely connected with each body influencing the other two. Therefore he cleaned and developed 
all three bodies for best results. 

“Philosophy” comes from the Greek words “philein” (meaning “to love”) and “sophia” (meaning 
“wisdom”). “Mathematics” comes from the Greek word “mathema” (meaning “science”), which origi- 
nates from “mathesis” (meaning “knowledge”). Therefore, in the true sense of the word, mathematics 
is the only science we have . 6 Pythagoras said: “Every man has been made by God in order to acquire 
knowledge and contemplate. ” So we should do mathematics in order to acquire knowledge and we 
should contemplate in order to transform our knowledge into wisdom. This makes mathematics the 
path to philosophy. 

Next about freedom. According to Pythagoras, the true world of our mind is non-material. However, 
because of our physical body our mind collects only material experiences. Therefore, our mind is “im- 
prisoned” by the material. Physically we never can be free, nor can we be free emotionally. But men- 
tally we can be free, if we break out of the chains of our mind’s material imprisonment. Mathematics 
can help with this, because the objects of mathematics are between material and immaterial. 

Look at the example of a point. As a mathematical object, a point is infinitely small. But there is noth- 
ing like this in the natural (material) world. However, there are objects (such as a point drawn with a 
pen, or an atom or subatomic particle) which come close to it for they are very, very small. Another 
example is a line, which, in mathematics, is infinitely thin and infinitely long. Also such a line does 
not exist in the material world, but there are objects (like a line drawn with a pencil and a ruler on a 
piece of paper) which comes somehow close to it. When we talk about a triangle in mathematics, we 
have something “in mind”, which can be considered an abstraction of triangular objects in natur or 
drawn on paper. Alternatively we can say that a triangle which we draw is a realization of a mathe- 
matical triangle. 

Immaterial 

ra 

1 I /Mathematics 


Ml 


Material 

Mathematical objects are like a rung on a ladder from the material world to the immaterial world. 
Therefore, doing mathematics helps our mind to shake off the confinements of the material experi- 
ences and raise above the material world. 


6 Natural sciences are sciences insofar as they use the methods of mathematics. 
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Now you probably wonder what this has to do with “numerics” and “symbolics”. 

Pythagoras said: “Number is the within of all things. ” The “things” are the objects in the natural (ma- 
terial) world. For Pythagoras the numbers are representations of these natural objects. This is why we 
call them “natural numbers”. The “natural numbers” are the most material of all mathematical objects 
insofar as they are very close to material things. (“5” is a mathematical object very close to the five 
fingers of a hand.) From the natural numbers we construct new objects which may be less material, 
such as the negative numbers. (“-5” is not so easily recognized in the material world.) 

Putting parts (remember: “part” = “number”) together as a new whole (remember: “whole” = “sym- 
bol”) is a basic technique in mathematics which leads to “less material” = “less numeric” = “more 
symbolic” = “more immaterial” objects. 


Immaterial 



Symbolics 

Numerics 


Material 


If we look at mathematics as (part of) a ladder from the material world to the immaterial world, “nu- 
merics” is the lower end (the first rung) on it, “symbolics” are the higher rungs. This means that the 
higher a person’s mathematical education, the higher on the ladder (s)he reaches. If you are high on 
the ladder there is the danger to lose contact to he material world. The picture of an absentminded 
mathematics professor appears ... However, the really great minds are those who are “tall” enough to 
stand with both feet on the material ground and reach with their hands high on the ladder. 

This picture is also supported by the following quote from Nikolaus of Cues Cusanus, a German theo- 
logian and humanist who lived 1401-1464: “If there is no other path to the divine then through sym- 
bols, we should use mathematical symbols for they possess indestructible certainty. Knowledge about 
the divine is out of reach for the mathematically illiterate. ” 

The following table summarizes our philosophical mediation: 


number 

symbol 

more material 

less material 

material 

immaterial 

finite 

infinite 

physical 

spiritual 


When “number is the within of all things ” (Pythagoras), then symbolics can be considered the within 
of mathematics. The Scottisc mathematician Eric Temple Bell (1883 -1960) said: “Any impatient 
student of mathematics or science or engineering who is irked by having algebraic symbolism thrust 
upon him should try to get along without it for a week. ” 
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We can’t do mathematics without symbolics, but also we shouldn’t do mathematics without numerics 
either for the numbers connect mathematics with the real world. 


Closing Remark 

I do hope that you found the contents of this presentation trivial, because, speaking with C F Gauss: 
“When a philosopher says something that is true then it is trivial. When he says something that is not 
trivial then it is false. ” 
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Appendix 1 


#1: nodi := 1 

intgr(x) := 

If INTEGER? 00 
x 

#2: If x > 0 

FLOOR (x) 

FLOOR 00 + 1 

nodi 

intgr(x-10 ) 

#3: r(x) := 

nodi 

10 



( ( 1 '\ 


#4: 

testl(n) := r r 

■r(n) 


L L n j 

j 

#5: 

testl(n) 



Appendix 2 


#1: nodi := 1 


intgrOO := 

If INTEGER?(x) 
x 

#2: If x > 0 

FLOOR (x) 

FLOOR (x) + 1 

nodi 

intgr(x-10 ) 

#3: r(x) := 

nodi 

10 


#4: test(x, y) := Lr(r(x) ) - r(r(y) ), r((r(x) + r(y)MKx) - r(y)))J 

2 2 

#5: testdiff (x , y) := r(r(x) ) - r(r(y) ) - r((r(x) + r(y))-(r(x) - r(y))) 

#6: testdiff (x, y) 
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Ein Brief von Mandala / A Letter from Mandala 


Sehr geehrter Herr Bohm 

Vielen bank fiir Ihr E- mail und Ihre Gluckwiinsche fur meine Aufnahme in der Open University. 
Auch Ihre Hinweise zu MathCad sind sehr niitzlich! Ich werde sie gerne beherzigen. Wie es 
scheint, ist diese Software integraler Bestandteil des Mathematik-Kurses der Undergraduate- 
Stufe der OU, so dass ich damit ein neues Programmiersystem kennen lerne. Umso besser! 

Meine Bewerbung bei der OU hat sich in der ersten Phase mithilfe von Korrespondenz ohne Kopie 
vollzogen, da aufgrund der OU- Satzung ich mir nicht im Entferntesten - das Aufnahme-Min- 
destalter betrtigt 18 Jahre- eine Chance ausgerechnet hatte, Student der OU zu werden. Das 
zweieinhalbstundige Kollogium, dass sich um meine Englisch - und Mathematikkenntnisse drehte 
(insbesondere um meine Fertigkeit, rekursive Programme mit dem TT- 92 und der Turtle- Geo- 
metry zu erstellen) wurde per E- mail anberaumt. Der leitende Mathematik- Tutor versicherte 
mir alsdann, dass mein Paper im DERIVE- Newsletter liber die Turtle- Geometrie ihn sehr bein- 
druckt und dieses auch wesentlich zu meiner Aufnahme in die OU beigetragen habe. 

Ich kann jedoch gerne bei ntichster Gelegenheit eine schriftliche Bestatigung des Vorgangs fiir 
Sie bei der OU an fordern. 

Auch ich wiinsche Ihnen einen scho'nen Sommer 
Mandala 


You might remember Mandala von Westenholz's contribution on TURTLE Graphics with the TI-92 
from DNL#57. Mandala writes in her mail that she applied for acceptance in Open University. She did 
not expect being accepted because minimum age must be 18. But she was invited for a 2 1/2 hour 
colloquium testing her English and Mathematics knowledge (recursive programming a.o.). The leading 
math-tutor told her that especially her paper published in the DNL had impressed him very much and 
that this paper and that this paper was responsible for a great deal that Mandala is now student of the 
Open University. Much luck, Mandala for your future and it would be great to hear from you. Thanks 
also to Bernhard Kutzler who forwarded Mandala's contribution and suggested publication. Josef 


Von Josef an Hans-Jiirgen Kayser 

At the occasion of revising and rewriting DNL#6 I wrote a letter to Hans-Jiirgen Kayser concerning 
his contribution to "Graphic Differentiation". In DOS-DERIVE times it was possible to slow down the 
plotting procedure by setting Plot Accuracy. So we could produce "mathematical movies". DERIVE 
and the PCs of our days are so fast, that it is not possible to reproduce this nice effect by means of 
1992. I didn't want to leave this "movie effect" and found a little trick to implement a "DELAY- 
function". 

Hans-Jiirgen Kayser answered very soon and thanked for "saving" the movie effect. He wrote that he 
also dealed with this effect by writing his Calculus Book (find a lot of H-J’s publications on the book- 
shelf). He found a similar way to resolve this problem, but uses now a mixture of his idea and my idea 
(using the RANDOM-function). He added some files to experience DELAY in Action. 


See one example on the next page. 
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You can find detailed instructions in "Derive im Analysisunterricht der Jahrgangsstufen 12 und 13", bk-teachware, SR-45 

#1: [CaseMode := Sensitive, In put Rode := Word] 

#2: [K(r , t) := [r.COS(t), r-SINCt)], n(r, t) := [r-t, r]] 

2 2 2 

*3: DELAY (dur) := VECTOR((x - RANDOM (1)) + (y - 200) = k_-0.02 , k_, 1, dur) 

#4: Zyk"loide_a(r , a, t) := [r-t - a-SIN(t), r - a-COS(t)] 

f 71 

#5: GRAFIK(r, a, d := 1) := VECTOR [K(r, T) + Pl(r, t) , [Pl(r, t) , Zyk1oide_a(r , a, t)], DELAY (d) ] , t, 0, 5 -jt, 

L ia 


f ” 

#6: ZYKLOIDECr, a) := VECTOR Zyk1oide_a(r, a, t) , t, 0, 5 -tt, 

L iao 

#7: GRAFIKQ, 1) 

*a: ZYKLOIDEC1, 1) 



The third parameter 3 delays the animation (default = 1) 

#10: GRAFIKQ, 2, 3) 

/II : ZYKLOIDEQ, 2) 



#12 : GRAFIKQ, 0 . 5 ) 

#13 : ZYKLOIDEQ, 0 . 5 ) 
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Walter Wegscheider 


Lieber Josef, 

ich bin gefragt worden, wie die Funktion der Spirale lautet, die 
DERIVE im Cover der CD verwendet. 1st anscheinend eine 
Art Mischung aus Torus und Spirale. WeiBt du, wie die Funk- 
tion dazu lautet? 

Walter was asked for the function to create the 
"DERIVE-SPIRAL" which is shown on the CD-Cover and 
DERIVE Manual 

Here is a collection of snails: 



#1: N0RMAL_VECT0R(v, t) := SIGN 


'd J2 


#2: B I NORMAL (v , t) := 5 IGn|cR05s| — v , - j v 

#3: snail (v , t, r, 4 .) := v + r- 1- (SIN C«t>3 ■ BINORMAL (v , t) + COSOtO ■ N0RMAL_VECT0R(v , t)) 

Dabei ist v die Parameterdarstellung einer Schraublinie mit dem Parameter t 
v is the parameter form of the helix. 

#4: snail ([S-COSOO , S-SINft), 0.9- t], t, 0.16, s) 

#5: snail C[5-C0S(t) , 5-SIN(t), 0.6- t], t, 0.1, s) 




Let's have an animated "spiral snail" with slider bars for a, b, c and k: 

r b-t b-t 1 

#6: spfr := |_a-e -COS(t) , a-e ■ SIN ( t]) , c-tj 

snai "l(spir, t, k, 5) 



# 7 : 
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Here is another mail from our colleague from Diisseldorf, Hans-Jiirgen Kayser. Once more 
I'd like to draw your attention on his inspiring publications (Information page). Josef 

Lieber Herr Bohm, 

anbei eine kleine Datei, die das Ziehen mit Zuriicklegen (vgl. DNL # 52, Seite 42) m. E. etwas einfa- 
cher lost. Ich hoffe, dass Sie sich hierfur interessieren. 

HJK 

Dear Mr. Bohm, 

I attach a short file, which in my opinion solves simulation of "Drawing without Repetitions" (DNL#52, 
page 42) a bit easier. I hope that you will find it interesting. 

Hans-Jiirgen Kayser 


Lotto / Lottery 

Simulation des Zahlen-Lottos "k aus n" 

Kurzbeschreibung / Short Description 


©HJK 6/2005 


Es ist kein Problem, Derive aus einer Menge von Zahlen (z.B: {1, 2, ..., 49} ) 6 Zufallszahlen auswah- 
len zu lassen. It is no problem for Derive to choose 6 random numbers from a set of numbers, eg from 
{1,2, ..., 49}. 


#1: AUSWAHL(n, k) := VECT0R(1 + RANDOM(n) , j, 1, k) 

#2: AUSWAHL(49, 6) 

#3: [17, 29, 12, 26, 46, 48] 

#4: [19, 26, 29, 31, 19, 29] 

Wie die zweite Serie zeigt, konnen bei den so definierten "Ziehungen" Wiederholungen auftreten. 

Die Funktion AUSWAHL simuliert das "Ziehen mit Zuriicklegen", da bei jeder neuen "Ziehung" wieder 
alle n Zahlen zur Auswahl stehen. 

Beim Zahlenlotto wird aber ohne Zuriicklegen gezogen. Dieses "Ziehen ohne Zuriicklegen" soli im 
Folgenden simuliert werden. 

As the second sequence of numbers shows repetitions of numbers can occur in the drawings defined 
by function #1 . In the following we will simulate "Drawings without Repetitions". 


Voriibung / Preparatory Exercise: 

#5: LOTTO 1 :z PR0G(P00L :z {1 49}, L00P(IF(DIM(P00L) = 43, 

RETURN {1 49} \ POOL), POOL := POOL \ {1 + RAND0M(49) })) 


Zwei Beispiel-Ziehungen / Two sample drawings: 

#6: L0TT01 = {12, 17, 26, 29, 46, 48} 

#7: L0TT01 = {22, 27, 29, 40, 41, 45} 
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Lotto "k aus n"\ 

LOTTOCn, k, POOL) :z 
Prog 

POOL :: {1 n} 

#8 : Loop 

If DIM (POOL) = n - k 

RETURN {1 n} \ POOL 

POOL := POOL \ {1 + RANDOM(n)} 


Hinweis: #8 ist folgendermaden einzugeben (vgl. # 5) / Hint: #8 must be authored as follows: 


LOTTOCn , k , POOL) : =PROG(POOL :={ 1 n} , LOOP(IF(DIM(POOL)=n-k , RETURN{ 1 n}\ 

POOL) , POOL : =POOL\{ 1+RANDOM (n) }) ) 


Beispiele / Samples: Lotto "6 aus 49", "3 aus 6" und "5 aus 12" 

#9: L0TT0(49, 6) = {3, 12, 15, 47, 48, 49} 

#10: LOTTOC6, 3) = {3, 4, 6} 

#11: LOTTO(12, 5) = {1, 2, 6, 11, 12} 


Eine Liste mit Tippreihen zum Lotto "6 aus 49" / A Sequence of drawings "6 out of 49": 


# 12 : 


#13: 


VECTOR( [ LOTTO(49 , 6)], k, 1, 10) 
' {12, 17, 26, 29, 46, 48} ' 

{22, 27, 29, 40, 41, 45} 

{1, 11, 26, 33, 46, 48} 

{4, 7, 22, 32, 36, 44} 

{1, 10, 14, 20, 22, 45} 

{11, 14, 22, 24, 28, 31} 

{30, 35, 37, 39, 42, 48} 

{1, 15, 20, 22, 24, 31} 

{4, 8, 13, 18, 30, 31} 

. {2, 9, 14, 21, 25, 49} . 


Alternative (aus: DERIVE-NEWSLETTER #52, p 42): 


LOTTERY(k, list, t_, newlist, counter) :: 
Prog 

newlist := [] 
counter :z 1 
Loop 

#14: If counter > k 

RETURN SORT(newlist) 
t_ :z RAND0M(DIM(1 ist)) + 1 
newlist := ADJ0IN(1 istj.t_, newlist) 
list := DELETEQist, t_) 
counter :+ 1 
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Beispiel-Ziehungen / Sample Drawings: 

#15: L0TTERY(6, [1 49]) = [11, 17, 26, 30, 46, 48] 

#16: VECTOR(LOTTERY(6, [1 49]), j, 1, 10) 


■ 11 

17 

26 

30 

46 

48 ' 

20 

26 

29 

40 

44 

45 

1 

26 

32 

40 

46 

48 

4 

8 

11 

33 

37 

45 

1 

10 

14 

20 

22 

46 

10 

14 

22 

24 

29 

32 

22 

31 

35 

38 

42 

48 

1 

19 

22 

24 

34 

42 

12 

14 

17 

23 

31 

33 

. 4 

8 

10 

15 

23 

30 . 


Simulation des Zahlen-Lottos "k aus n”l Ende 

I'd like to add a tip: Include simplifying RANDOM(O), either explicitly as first command or built in 
in the frame of the program (eg. dummy :=RANDOM(0)), to obtain different simultations at each 
new session. Josef 


Bill Martin bill_e_martin0HOTMAIL.COM 

Apologies if this is the wrong forum in which to pose this question. 

I use Derive 5 and sometimes need the Kronecker product of two matrices. Matlab has a "kron" func- 
tion but I believe Derive does not possess an equivalent. 

The Kronecker product can be formed in Derive using a longwinded manual method involving the 
commands “ELEMENT”, “APPENDCOLUMNS” and “APPEND”. 

I would much prefer to have a simple command function but do not possess the requisite skills to pro- 
gram in Derive. 

I would be most grateful if anyone could offer assistance, suggestions and, best of all, the Derive code 
for a Kronecker product function. 

Thank you. 

Bill, 

I don’t remember if Derive 5 has the MAP_LIST function, but it if does you could use something like 
this. 

Kronecker Product (a, b) := APPEND(MAP_LIST(APPEND_COLUMNS(x_) , x_, MAP_LIST 

(MAP_LIST(v_-b, v , w ) , w_, a))) 


Don Phillips 
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Titbits from Algebra and Number Theory (31) 

by Johann Wiesenbauer, Vienna 

Have you ever used sets and set-theoretic operations in your programs? If your 
answer is no, then you may have missed something. In fact, they come in handy 
on many occasions, where most people don't think of sets in the first place. 

Let's start with a very simple example of this sort. Imagine 30 hunters aiming at 
10 ducks in front of them and shooting at exactly the same time. We assume 
here that each hunter selects the target independently of the others and that 
there is no miss. What do you think is the chance that all 10 ducks were hit af- 
ter that shooting? Would you dare to bet on this event? 

Well, I'm sure most mathematicians will immediately translate this question into 
the following one: What is the probability that a random mapping from a 30- 
element set into a 10-element set is surjective? Now some will try to answer 
this question in a purely theoretical way, which is not too difficult, namely 

2-30 

APPROX(5TIRLING2(30 , 10).10!-10 ) = 62.91371892 

Here, STIRLING2(30,10) is the number of ways to decompose a 30-element set 
into exactly 10 nonempty subsets. These correspond to partitions induced on the 
30-element set by any surjective mapping onto the 10-element set, where ele- 
ments belong to the same class if and only if they have the same image under 
the mapping at issue. For each of these partitions you can select images of its 
members in 10! ways, hence the factor 10! above. Finally, one must divide 
STIRLING2(30, 1 0) 10! by the number of all mappings, i.e. 10 30 , and multiply it by 
10 2 to get the percentage. 

Many Derivers will prefer an "experimental" approach though, which is sketched 
below. First of all, we need a routine that determines whether a mapping from an 
n-element set into a k-element set is surjective or not. This is where our sets 
come into play, as you can see in the following implementation. 

sur ject"ive?(n , k, s_ := {}) := 

Loop 

If n = 0 

RETURN DIM(s_) = k 
s_ s_ u { RANDOM (k)} 
n 1 

Here and in the following, we make use of the basic fact about sets that every 
element in a set occurs only once and adding an element already contained in the 
set doesn't change it. 

The next routine performs a number of tests w.r.t. surjectivity (hence its uni- 
maginative name!) and after that computes the proportion of surjective map- 
pings as a percentage. 
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testOn, n, k, m_ := 0, t_ := 0) := 

Loop 

If m_ = m 

RETURN 100 ■ t_/m 
If sur ject"ive?(n , k) 
t_ :+ 1 
m_ : + 1 

APPR0X(test(10000, 30, 10)) =62.86 

In particular, we can see a good coincidence with the theoretical result and they 
both show that is worthwhile to bet on all ducks getting hit. Is this what you 
guessed? 

Our next example is from coding theory. A binary linear code C could be given by 
a matrix G with entries in Z2 along with the rule that a vector c is in C if and 
only c can be written as a sum of certain row vectors of G. For example, if 

' 1 0 0 0 0 1 

G:z 010100 

.001010 

then c = (1,1, 0,1, 0,1) belongs to the code, because it is the sum of the first two 
row vectors. How could a program look like that computes the whole code C? 

Would you believe that this small problem usually gives most Derivers quite a 
headache, even if I give as a hint that set-theoretic functions should be used? 
Ok, I don't want to put you on tenterhooks no more, here is the surprisingly 
short solution using Cartesian products. 

3 

C :z VECTOR(MOD(v-G, 2), v, {0, 1} ) 

' 0 0 0 0 0 0 
0 0 10 10 
0 10 10 0 
0 11110 

C := 

1 0 0 0 0 1 
10 10 11 
110 10 1 
111111 

If you prefer the representation of C as a set rather than a list you should use 

3 

C := MAP_LI ST (M0D(v ■ G , 2), v, {0, 1} ) 
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C := {[0, 0, 0, 0, 0, 0], [0, 0, 1, 0, 1, 0], [0, 1, 0, 1, 0, 0], [0, 1, 1, 1, 1, 0], [1, 0, 0, 0, 0, 1], [1, 0, 1, 
0 , 1 , 1 ], [ 1 , 1 , 0 , 1 , 0 , 1 ], [ 1 , 1 , 1 , 1 , 1 , 1 ]} 

After these warm-up exercises, it's high time we tackled a more ambitious pro- 
ject dealing with sets. I'm speaking of a routine that helps you a lot when solving 
Sudokus, the well-known Japanese cult puzzles. For members of the supposedly 
minority, who never heard of it, here a short explanation. Given is a 9x9-matrix 
like the matrix A below with some nonzero entries in the range 1-9. 

000007500 

040006300 

002300000 

013000000 

530000027 

000000960 

000004100 

007900030 

006200000 

For the solution of a sudoku it is required to replace all zeros, which are actually 
empty entries in a real sudoku, again by numbers in the range 1-9, such that 
every row, every column and nine 3x3-submatrices contain all numbers 1-9. As 
for the nine 3x3-submatrices, you will see them in the solution below, which was 
obtained by the program on the next page and the input below. 

A := sudokuQ 




19 6 


7 4 3 
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sudoku (n := m, a_, b_, c_ := 1, i_ := 1, j_ := 1, k_ := 1, 1_ := 1, s_, t_, u_) := 

Loop 

If DIM (A) = 9 
Prog 

A :+ VECTOR (VECTOR (123456789 ■ 0 a A| j , j , 9) , i , 9) 

A := VECTOR (VECTOR (VECTOR (VECTOR (Aj,(i + k)±(j + 1), 1, 3), k, 3), j, 0, 6, 3), 
If n = 0 v c_ = 0 
RETURN A 
n 1 
c_ := 0 
d_ := 1 
Loop 

a_ := ABS(Aj, i_ 4 , j_j,k_j,l_) 

If a_ = 0 

RETURN ''Computation has been stopped due to zero entries! 11 
s_ := MAP_LI ST (a_i.m_ , m_, {1 DIM(a_)}) 

t_ := TERH5a(x A AE5CAn_4(m_4l)|k_iCm_|2)) 1 m_, {1, 2, B} A 2 \ {[j_, !_]})) 
u_ := SUBST (t_, x, 1) 
t_ := x-SCLNCO, x) 

Loop 

Ifu_= [] exit 
If DIH(FIRSTCO) = FIRST (u_) 

s_ := s_ \ MAP_LIST ((FIRST m_, {1 DIM(FIRST(0)D 

u_ := REST (u_) 
t_ := REST CO 

t_ := TERHSCICx A ABSCAUm_tl)lJ_lCm_|2)|l_) 1 m_, {1, 2, B} A 2 \ {[i_, k_]})) 
u_ := SUBST (t_, x, 1) 
t_ := X-3CLNCO, x) 


Loop 

If u_= [] ex it 
If DIHCFIRSTCO) = FIRSTCO 

s_ := s_ \ HAP_LI ST C (FIRST Ct_)) j.m_, m_, {1 DIH(FIRSTCO) }) 

u := REST (O 

t_ := REST (t_) 

t_ := TERHS(I(x A ABS(A|i_4.j_j.(m_4.1)|(m_j.2)) l m_, {1, 2, 3} A 2 \ { [ k_ , !_]})) 

u := SUBST (t_, x, 1) 

t_ := x-5(LN(t_), x) 

Loop 

If u_ = [] exit 

If DIM (FIRST (t_)) = FIRST (u_) 

s := s \ MAP_LI ST ((FIRST (t_) )|m_ l m_, {1 DIM(FIRST(t_)) » 

u := REST (u_) 

t_ := REST (t_) 
t_ := 0 
Loop 

If s_ = {} ex it 

t_ := 10- 1_ + FIRST (s_) 

s := REST (s_) 

b_ := IF(a_ = t_) 
c : + 0 A b_ 

A.j. i_.|. i—i-k—i. 1_ := (2 ■ b_ - 1) ■ t_ 

1_ :+ 1 
If 1_ > 3 

[1_ := 1, k_ :+ 1] 

If k_ > 3 

[k_ := 1, j_ :+ 1] 

If ]_ > 3 

[j_ := 1, i- :+ 1] 

If i_ > 3 

[i_ 1, exit] 


6, 3) 


I don't want to explain this program in detail, but only point out some important 
points. First of all, the first parameter is sort of a switch, by which you can con- 
trol how many steps are performed before the program stops. (Note that, if 
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there are no changes anymore as to the entries of the matrix, it will stop any- 
way!) In particular, you can set up a step-by-step mode by choosing n=l. In this 
case, changed entries w.r.t to the last matrix will be indicated by a minus sign. 

For example, setting n:=l we get for the same matrix A as above 

A := sudoku(l) 



For sudoku solvers, who want rather a "hint" than the whole solution, this could 
also be what they want. Furthermore, you can observe, as it were, how the pro- 
gram “thinks", which might be interesting as well. 

Let's have a look at the very first entry -13689. It says that there was a 
change, namely only the digits 1,3, 6, 8, 9 are possible in this place instead of all 
digits before. Scanning the first row of the original matrix A, you see that 7 and 
5 are excluded, indeed. As for the first column, there is nothing else we can 
conclude from it. Finally, the first 3x3-submatrix says that 2 and 4 are also im- 
possible in top left corner, leaving only the possibilities 1,2, 6 ,8, 9, indeed. 

Now let's focus on the two entries -35 at the end of the fourth row. What do 
they tell us? Of course, experienced sudoku solvers know the answer: The digits 
3 and 5 are excluded in the rest of the row and also in the rest of the 3x3- 
submatrix. The same is true if a number with 3 digits occurs three times in a 
row, column or 3x3-submatrix and so on. 

Actually, this is the only strategy the program pursues. Hence, it could be that 
the program stops although it has not yet arrived at a solution. In this case it's 
up to you to edit the output matrix in some places (e.g. by trying out one of two 
possibilites) and the program will continue from there. Nice, isn't it? I hope you 
enjoy the program as much as I do. 
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As you can imagine, I immediately tried this tool - I am not a SODUKOer, but my wife Noor is and I 
found several Sodukos in the newespapers. Easy and Medium was no problem for the program, but 


EASY 


MEDIUM 


'DIABOLICAL" 



9 

6 


5 


3 

T 


4 



8 


3 



9 

1 




7 




4 


5 


6 


8 


i 


8 


2 




4 


6 


4 


9 


1 


3 


9 




1 




2 

5 



7 


2 



3 


3 

8 


9 


7 

5 





4 




7 





7 

5 

1 



1 



4 

6 

3 


2 


5 

2 




1 

3 


9 

7 




8 

4 

~ 

~r 

~3~ 

~ 


~ 

~e~ 

2 





~ 9 ~ 



3 




~2~ 


~e~ 





~8~ 







*2: A :z 


#4: A :z 


096050320 
400803009 
100070004 

0 5 0 6 0 8010 #3 : sudoku O = 

802000406 

040901030 
900010002 
500702003 
03 80907 5 oj 3 S3 sec 

004000700 
000751000 
100463002 
052000130 

#5 : sudoku 

097000840 
013000620 
500198003 
000246000 
0080004 00. 4 03 sec 

[ 1478 12347 5 

4678 3467 3478 

L 14678 9 13478 


#7: sudoku () = 


#5: sudoku () = 


7 9 6 

4 2 5 

.18 3 

3 5 9 

8 12 
.647 
'974 

5 6 1 

.238 

3 6 4 

2 8 9 

17 5 

8 5 2 

6 9 7 

4 13 

5 4 6 

7 3 1 

9 2 8. 




5 




6 






1 


2 





9 






2 


3 



5 


6 



1 

9 








7 

~2 

~ 


~ 3 ~ 













1 





~T~ 


~Y~ 






~ 6 ~ 




IT 




15 4 

8 6 3 

.279 

6 4 8 

5 3 7 

.921 
3 15 

7 8 2 

.496 

8 2 9 

7 5 1 

4 6 3 

6 7 4 

3 12 

9 8 5 

19 8 

2 4 6 

5 3 7 . 


3 2 8 

17 9 

.564 
2 17 

4 9 6 

.835 

6 8 2 

9 4 3 

.751 

7 5 1 

3 6 4 

9 8 2. 

13 9 

8 4 5 

6 2 7. 

2 7 3' 

5 9 8 

4 16. 


3 

47 

478 

9 

1456 

148 

2 

567 

78 


457 8 23479 

157 1257 1279 

1457 123457 6 


789 34789 4789 

1 3456789 2 

678 345678 4578 

5 24789 6 

28 248 48 

. 3 789 1 . 

2679 25679 579 

4 1256789 3 

279 12579 579 


6 34789 389 

34579 345789 3589 

. 13457 2 358 

29 89 1 

235 3568 7 

. 59 5689 4 . 

234579 1 23569 

2579 5679 2569 

8 34579 2359 


1 .78 sec 
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The ..DIABOLIC" was not so easy to solve. Johann sent a mail with some hints referring to the 
last paragraph of his contribution. 

Program terminates after 2 iterations. We proceed by systematically distinguishing possible 
cases: investigate systematically possible cases. So on position (6.3) only 7 or 8 are possible. We 
issue 7 and simplify A:=sudoku() - without pressing RETURN, but using the = command, either by 
pressing the —Button left of the entry-line. 7 on this position leads to a contradiction. 


*7: 


1478 12547 5 

4678 5467 5478 

14678 9 15478 




5 

47 

478 


sudoku () = 


9 

1456 

148 




. 2 

567 

m . 




' 457 

8 

25479 



157 

1257 

1279 



. 1457 

125457 6 


* 8 : 


A := 


1478 

12547 

5 

4678 

5467 

5478 

. 14678 

9 

15478 


5 

47 

478 



9 

1456 

148 



. 2 

567 

7 . 


457 

8 

25479 

157 

1257 

1279 

. 1457 

125457 

6 


789 54789 4789 

1 5456789 2 

678 545678 4578 

5 24789 6 

28 248 48 

5 789 1 

2679 25679 579 

4 1256789 5 

279 12579 579 

789 54789 4789 

1 5456789 2 

678 545678 4578 

5 24789 6 

28 248 48 

5 789 1 

2679 25679 579 

4 1256789 5 

279 12579 579 J 


6 54789 589 

54579 545789 5589 

15457 2 558 

29 89 1 

255 5568 7 

59 5689 4 

254579 1 25569 

2579 5679 2569 

8 54579 2559 

6 54789 589 

54579 545789 5589 

15457 2 558 

29 89 1 

255 5568 7 

59 5689 4 

254579 1 25569 

2579 5679 2569 

8 54579 2559 


#9: sudoku () = Computation has been stopped due to zero entries! 


Taking 8 on this position leads to the next step, eg, checking 1 and 2 on position (1,2): 


18 

EE 

5 


789 

54789 

4789 


6 

54789 

589 

68 

547 

547 


1 

5456789 

2 


54579 

545789 

5589 

. 168 

9 

547 . 


. 678 

545678 

4578 . 


. 15457 

2 

558 . 


5 47 47 


5 289 6 


29 89 1 

9 56 1 


28 248 48 


255 5568 7 

. 2 56 8 . 


. 5 79 1 . 


. 59 569 4 . 


457 

8 

25479 


2679 

25679 

579 


254579 

1 

25569 

157 

127 

279 


4 

1256789 

5 


2579 

5679 

2569 

. 1457 

12547 

6 


. 279 

12579 

579 . 


8 

54579 

2559 . 


2 turns out to be the right choice. 

So you might go on. You are invited to download the solution. It is among the Derive files. 







